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I. INTRODUCTION 



The seminal paper by Bender and Boettcher in 1998 has lead to an alternative formulation 
of quantum mechanics (QM): non-Hermitian PT-symmetric QM [1]. Although the Hamil- 
tonians involved appear to be not Hermitian {H ^ H'^), they yield only real spectra and the 
time evolution is unitary when the VT symmetry is not broken 

Since the introduction of PT-symmetric QM, there has been a debate about whether 
PT-symmetric QM is more general than the conventional QM. In the first paper about 
VT symmetry, Bender and Boettcher stressed that VT symmetry is a weaker condition 
than Hermiticity and a PT-symmetric theory can be considered as a complex extension of 
a conventional Hermitian theory In 2002, Bender et al. constructed a positive-definite 
norm of PT-symmetric QM using the newly discovered C operator p]. The claim that 
VT symmetry is more general than Hermiticity was restated in the paper's title: "Complex 
Extension of Quantum Mechanics. " Interestingly, in the same year, Mostafazadeh pointed 
out that all pseudo-Hermitian Hamiltonians are Hermitian (self-adjoint) with respect to a 
positive-semidefinite inner product [s*]. As such a PT-symmetric Hamiltonian can be con- 
sidered as a special case of a pseudo-Hermitian Hamiltonian with a real spectrum. In 2003, 
Mostafazadeh showed that diagonalizable pseudo-Hermitian Hamiltonians are extensions of 
VT- or CPT-symmetric Hamiltonians . A less general result was obtained independently 
by Bender et al. in Ref. 0]- But they took a different point of view which can be seen from 
the title of the paper: "All Hermitian Hamiltonians Have Parity. " 

The study of finite-dimensional matrix Hamiltonians may shed some light on the debate. 
The original 2x2 "PT-symmetric matrix Hamiltonian introduced in Ref. [2| has three real 
parameters. In 2003, Bender et al. extended it to a four-parameter class and to higher 
dimensions jl] . From their results, they identified PT-symmetric Hamiltonians and Hermi- 
tian Hamiltonians as two distinct extensions of real-symmetric Hamiltonians. For matrix 
dimension higher than two, PT-symmetric Hamiltonians have less real parameters than 
Hermitian Hamiltonians in the same dimension. This view was challenged immediately by 
Mostafazadeh in a paper titled "Exact VT -Symmetry Is Equivalent to Hermiticity" [i9|. In 
this paper, the author constructed a 2 x 2 PT-symmetric Hamiltonian with five real pa- 
rameters. In 2007, Mostafazadeh and Ozgelik explicitly constructed the most general 2x2 
quasi- Hermitian Hamiltonian with six real parameters [lo| . 

In this paper, we are trying to settle the debate by constructing the most general PT- 
symmetric matrix Hamiltonians. Part of the reason for the debate is due to the different 
usage of the terminology. To avoid further confusion, we restrict the term "Hermitian 
conjugate" or "f" in the Dirac sense: complex conjugate and transpose. In the present 
paper, we formulate PT-symmetric QM slightly different from the literature in two aspects. 
First, we define the time reversal operator as the Dirac conjugation rather than just the 
complex conjugation used in Ref. [H, 0, [8|. Second, we do not introduce any biorthonormal 
basis as in Ref. [sl, 0] . We link the CPT- inner product to the general inner product in a 
linear vector space with a weight function. The latter is the standard notation which can 
be found in modern quantum mechanics textbooks such as Ref. 

With this formulation, we solve for the general P operator and present the solutions 
explicitly in the case of 2 x 2 and 3x3. Using the general P operators, we construct 
the general PT-symmetric matrix Hamiltonians. We confirm that the most general 2x2 
PT-symmetric Hamiltonian has six real parameters as shown in Ref. |tlO|]- We find that the 
general 3x3 PT-symmetric Hamiltonian has thirteen real parameters. In the case of 2 x 2, 
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the general PT-symmetric Hamiltonian represents the most general matrix Hamiltonian 
with a real spectrum. Interestingly, this is not true in the case of 3 x 3. In both cases, 
we show clearly that all Hermitian Hamiltonians are just special cases of PT-symmetric 
Hamiltonians. From these finite dimensional results, we conjecture that VT symmetry is a 
generalization of Hermiticity in general. 

The paper is organized as follows. In Sec. [Tll we give our formulation of "PT-symmetric 
QM. In Sec. [ml we illustrate our ideas in the case of 2 x 2. All relevant operators or matrices 
are calculated explicitly. Special cases discussed in the literature are analyzed. In Sec. [TVl 
we construct the V operator and the PT-symmetric Hamiltonian in the case of 3 x 3. Finally, 
in Sec. |Vl we give some concluding remarks. 

II. PT-SYMMETRIC QUANTUM MECHANICS 

In this section, we formulate PT-symmetric QM. We start with a brief summary on the 



inner product in QM. We adopt the notation used in Ref. [11 



A. Inner Product 

In quantum mechanics, the Hilbert space can be considered as a linear vector space associ- 
ated with an inner product. The inner product between two quantum states, denoted {■■,■), 



must satisfy the following conditions |11| : 

1. (■?/', 0) is a complex number, 

2. ('?/', 0) = {(p.Tp)*, where * denotes complex conjugate, 

3. {ipi Ci4>i + C202) = Ci^ip, + C2{ip, (P2), where Ci and C2 are complex numbers, 

4. {4>, <f)) > 0, with equality holding if and only if = 0. 
In general, we may define the inner product as 

(^11^10), (2.1) 

where W is the weight function and the bra state is defined as the Hermitian conjugate of 
the ket state, (-1 = |-)^ [ll|]. 

From Eq. fl2.1l) and the first three properties of the inner product, it can be easily shown 
that W must be a Hermitian operator: W = W'^ . From the fourth property of the inner 
product, W has to be positive definite. That is, all the eigenvalues of W are positive. 

In the Hilbert space defined above, a self-adjoint operator, such as the Hamiltonian if, 
satisfies 

(^,/70) = (if^,0) (2.2) 
for arbitrary states and ip. All eigenvalues of a self-adjoint operator are real. And the 



eigenstates corresponding to different eigenvalues are orthogonal [11 . 

Plugging the self-adjoint condition (12. 2p into the definition of the inner product in (12.11) . 
we obtain 

WH = H^W. (2.3) 

We may consider this equation as the definition of the weight function W for a given Hamil- 
tonian H. Thus, the inner product is dynamic (Hamiltonian dependent) in general. 
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In conventional QM, the weight function can be chosen as the identity operator. In 
this case, the self-adjoint condition in (12 .Sp reduces to the Hermiticity condition: H = 
Since the identity operator is independent on Hamiltonians, the inner product is no longer 
dynamic. 

B. VT Symmetry 

We define the time reversal operator, T as Dirac conjugate. That is, for an operator A, 

rAr=AK (2.4) 

It follows that = 1, here U is the identity matrix. Note that our definition of the T 
operator differs from Refs.[l|, 0]. This definition allows us to have a more general parity 
operator. 

For the parity operator V, we demand it to commute with the time reversal operator and 
to be an involution. That is, 

1. [P, 7^ = 0, or equivalently, V = V'^ , 

2. p2 = 1 

There are obviously two trivial solutions to these constraints, 

Vo = ±1. (2.5) 

We will discuss the non-trivial 2x2 solutions of V in Sec. IIIII and 3 x 3 in Sec. IIVI 
A PT-symmetric Hamiltonian commutes with the combination operator VT: 

[H,Vr]=0 ^ VH^V = H. (2.6) 

Using the trivial solutions of Vq in (12. 5p . we get that a Po'7~-symmetric Hamiltonian is 
Hermitian 

[H,VoT\ = Q ^ H = Hl (2.7) 

In this sense, a Hermitian Hamiltonian is Po7~-symmetric with Vo to be the plus or minus 
identity matrix. 

In general, the eigenstates of a PT-symmetric Hamiltonian with different eigenvalues are 
not orthogonal with respect to the Dirac inner product. To solve this problem, one may 
define a VT inner product as 

(^,0)pr=(^l^l0). (2.8) 

From the commutation relation between H and VT, it can be shown that the eigenstates 
of H with different eigenvalues are orthogonal with respect to the VT inner product [H, [3|. 

However, the norm with respect to the VT inner product is not positive definite. This 
is simply because V has negative eigenvalues. One has to normalize the eigenstates to ±1. 
To overcome this difficulty, one needs to find a positive-definite norm by introducing the C 
operator [2]. Here, we define the C operator as 

C = J2\E^){E,\V, (2.9) 

i 

where \Ei) are the eigenstates of H with the PT-norm +1 or —1. 
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Our definition of tlie C operator lias tlie same properties as tfie one constructed in Ref. jl]. 
From tlie orthogonality and the (non-positive-definite) normalization of the Vl~ inner prod- 
uct, it can be shown that \Ei) are eigenstates of C with eigenvalue equal to the PT-norm: 

C\Ei) = {Ei\V\Ei)\Ei). (2.10) 

Thus, C commutes with H. The C operator also commutes with VT: [C, VT] = 0. This 
fact can be verified by using VC^V = C. Because the eigenvalues of are all unity, the C 
operator is an involution just like the V operator. 

Note that we do not require \Ei) to be simultaneous eigenstates of VT and of H. In fact, 
this can only be achieved in the special cases of symmetric matrices, such as Hamiltonians 
constructed in Refs. [1, [1] . We will discuss more details about these two examples in SecJIIIi 

Equipped with the C operator, we are ready to construct an inner product with the 
positive-definite norm. We define the CVT inner product as 

{ij,4>)cPT={i^\VC\(P). (2.11) 

Comparing to the general inner product in fl2.ll) . we recognize the weight function for the 
CVT-inner product is 

W = VC. (2.12) 

Since the weight function is positive-definite, we may find its square root W = iff, where 
7] is Hermitian, t]'^ = rj. Using the operator rj, we may define a Hermitian Hamiltonian h, 
which has same spectrum as a PT-symmetric Hamiltonian H j^, 

h = 'r]Hrj^^, where h = hK (2.13) 

We would like to emphasize that the transformation rj is Hermitian rather than unitary, 
hence the above relation is not unitary equivalence in the usual sense. 



III. 2x2 CASE 



In this section, we illustrate the ideas in the previous section by using 2x2 matrices. We 
find that any 2x2 Hermitian matrix is a special case of the general "PT-symmetric matrices. 
Furthermore, it is found that the most general matrix with a real spectrum must coincide 
with the general "PT-symmetric matrix we constructed. 

Since the V operator is an involution, it is a square root of the identity matrix. In 2 x 2 
matrices, other than the trivial roots in (12.51) . there is a non-trivial root with the form 

/ cos^ sin^e-i'^\ , . 

' l^sin^e'"^ -cos^ J' ^ 

where 6 and if are two real parameters. In terms of the Pauli matrices, the V operator can 
be written as P = n'' ■ cr, where n'" = (sin 6 cos if, sin 6 sin ip, cos 6) is a unit vector. 
To find the general PT-symmetric Hamiltonian, we use the following ansatz: 

H = et + a - cr, (3.2) 

where e and a = {ax,ay,az) are complex numbers. Plugging the above ansaiz into (12.61) . 
we get the equations satisfied by e and a: 

e = e*, a + a* = 2 (a* ■ n") n''. (3.3) 
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The first equation simply says that e is reaL The second equation can be written as 

^Mfc^a, = a*, with Mki = -5ki + 2nln\. (3.4) 

i 

In this form, the searching for a PT-symmetric Hamiltonian becomes an eigenvalue problem. 
If we separate the real part and the imaginary part as o: = A + i B, then A and B can be 
considered as the eigenvectors of the matrix M with eigenvalues +1 and —1, respectively. 
The matrix M always has one eigenvector with eigenvalue +1, and two eigenvectors with 
eigenvalue —1. It is easy to show that the eigenvector with positive eigenvalue is parallel 
to n'' and the eigenvectors with negative eigenvalue are perpendicular to n'^. Thus, the 
solutions to the eigenvalue problem are 

A = 7n^ B = /in'' + z/n'^, (3.5) 

where 7, /x, and v are real parameters and = (cos^ cosy9, cos 6* siny9, — sin6') and n*^ = 
(— sin cos 0) are two unit vectors. 

Thereafter, plugging ol = 'yif + i/in^ + iz/n'^ into the ansatz in (13. 2p . we get the general 
PT-symmetric 2x2 Hamiltonian matrix. 



H = el+ (711'^ + i/xn^ + iz/n^) ■ cr 

e + 7 cos 6 — ifi sin 6 (7 sin 6 + ifi cos 6 + v)e~^'^ 
(7 sin 6 + in cos 6 — z/)e"^ e — 7 cos 6 + 1^ sin 6 



(3.6) 



This Hamiltonian has six real parameters: e, 7, /i, z/, 9, and ip. All 2x2 Hermitian matrices 
can be recovered as special cases with fi = u = 0. In other words, the VT symmetry is a 
generalization of Hermiticity. 

The eigenvalues of the Hamiltonian H are 

E± = e± v/72 - /i2 - z/2. (3.7) 

For eigenvalues to be real, it requires 7^ > yU^ + z/^. For simplicity, we only consider the 
non-degenerated case with 7^ > /i^ + z/^ in this paper. The corresponding eigenstates are 



u 



3i«o-iv /I + £ sin^ ± ^^72 - _ ^2 cos 9 
V2 I e^-i /i _ £sin^T iv^72 -/i2 - u^cos9 



where we have defined three angles and a normalization constant as 

kq = arg (7 sin 6* + z/ + ifi cos 9) , 

K± = arg ^—7 cos 9 ± A/72 — /i2 — i/2 _|_ gj^ 5) j ^ 

''^J^^2 2- (3-9) 

y 72 — — z/2 

The normalization is chosen such that 

(E±|P|E±) = ±sign(7). (3.10) 



6 



The six-parameter-class of matrix in (13. 6 p with the condition 7^ > /i^ + coincides with 
the most general 2x2 matrix with only real eigenvalues. Qualitatively, this can be seen 
from parameter counting. A general complex 2x2 matrix has eight real parameters. The 
reallity of all eigenvalues puts two constraints on the matrix. Therefore, the most general 
2x2 matrix with only real eigenvalues should consist of six real parameters. 

This coincidence can also be proved rigorously. By direct computation, the eigenvalues 
of an arbitrary 2x2 matrix with the form of the ansatz in (13. 2p are E± = e ± ^/a~~a.. 
Imposing reality condition on the eigenvalues leads to 

e = e\ a-a>0. (3.11) 

The first condition is the same constraint on the parameter e as in the PT-symmetric 
Hamiltonian in (13. 6p . The second condition in (13. lip implies that the real part and the 
imaginary part of a are perpendicular to each other, A-B = 0, and that the real part vector 
is not shorter than the imaginary part vector, A ■ A > B • B. Without loss of generality, 
we may parametrize the real part vector as A = 711^'. Then the above conditions lead to a 
unique solution for the imaginary part vector which can be parametrized as B = yun^ + z/n'''. 
Since A is not shorter than B, we have 7^ > /x^ + z/^. Clearly then, the PT-symmetric 
Hamiltonian in (13. 6p represents the most general 2x2 matrices with only real eigenvalues. 

For the PT-symmetric Hamiltonian H in (13. 6p . the C operator can be calculated directly 
from its definition in (12. 9p . It is also straightforward to construct it from the relation in 
(I2.10p . Either way, the C operator is found to be 

C = -cx-(T. (3.12) 

7 

The form of the C operator is not a surprise because it is defined as an involution. Therefore, 
it must be a square root of the identity matrix just like the V operator. This fact can be 
easily verified by observing that {u/'~f)ct is a unit vector. 

Because the C operator has eigenvalues C\E±) = ±sign(7)|i?-|-), the eigenstates in (13. 8p 
are normalized to unity with respect to the CPT- inner product. Thus, we have a set of 
orthonormal eigenstates, 

{Ei\VC\E,) = z,j = ±. (3.13) 

The weight function has the form 

W = VC = u{l+f3-fT), where /3 = -n^ - -n^. (3.14) 

7 7 

Note that /3 is a unit vector and it is perpendicular to a. Interestingly, the square root of 
W has the form 

V2(m± 1) 

There are two solutions for rj, which is not because of the arbitrary overall sign in the 
square root. Rather, it is corresponding to two choices of mapping the eigenstates during 
the similarity transformation in (I2.13p . In the Hermitian limit, W ^ 1, rj is once again a 
square root of the identity matrix. The "+" sign in (I3.15P has the limit rj^ 1 and the 
"— " sign has the limit rj- ^ f3 ■ cr. 

Using the operator r], we find the Hermitian equivalence of PT-symmetric Hamiltonian 
H in dSSD, 

h = et±-n'' ■ a. (3.16) 
u 
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Now let us consider some special cases. We show that both Hermitian Hamiltonians and 
several PT-symmetric 2x2 Hamiltonians studied in the literature can be reduced from our 
general PT-symmetric Hamiltonian in (13. 6p . 



A. Special Case 1: Hermiticity 

If we set iJ, = u = 0, H becomes Hermitian, H = W, 

Tj ir , r / e + 7 COS 6^ 7 sin 6^ e""^ \ lo 

i^Hern^itian = 5l + 7n ■ ^ = ^ ^ ^.^ ^ . - 7 COS ^ J " ^^.17) 

This matrix Hamiltonian has four real parameters, and it includes all 2 x 2 Hermitian 
matrices. 

In this case, the weight function reduces to the identity matrix and the C operator coin- 
cides with the parity operator: 

W^Hermitian = ^Hermitian = V. (3.18) 

All these observations are consistent with the conventional QM. We may say that the 
Hermitian Hamiltonian is a special case of the PT-symmetric Hamiltonian with C = V. Or, 
equivalently, VT symmetry is a generalization of Hermiticity. 



B. Special Case 2: Bender-Brody- Jones Hamiltonian 

In Ref. [2], Bender et al. studied a "PT-symmetric Hamiltonian with real and symmetric 
off-diagonal matrix elements. Their choice of the parity operator is Pbbj = <^x- This case 
can be reduced from our general case by setting z/ = = and = -n /2. In particular, we 
have 

-=-(^v%;i.)^ ^e.^/g:(f ,i). (3.19) 



These expressions are equivalent to those in Ref. [2i by mapping our parameters e, /x, and 
7 to r cos ^, — r sin 9, and s therein. 

C. Special case 3: Bender-Meisinger-Wang Hamiltonian 



In Ref. [81, Bender et al. generalized the Hamiltonian matrix in Ref. [2j] by choosing a one- 
parameter class of parity operator, 

Pbmw = ""211 -cole ) • ^^-2°) 

This can be recovered by setting (f = in (13. ip . If we consider only symmetric Hamiltonian 
as in Ref. [sl], we can further set z/ = in (13. 6p . With this choice of parameters, we have 



BMW 



£ + 7 COS 6 — ifi sin 6 7 sin 6 + cos 6 
7sin6' + iyUcos6' £ — 7 cos + i/i sin 6' 




, , , cos6' — i- sin6' sin^ + i-cos^^ , 
'^BMW = \ I - ^ 1^ sin ^ + i^ cos 9 -cose +V smO ' ' 
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Once again, these are the same formulas as in Ref. [8| by properly mapping the parameters. 

The symmetric Hamiltonian, H in (13. 6 p with i/ = y9 = 0, has additional properties. In 
this case, the eigenstates of H are also the eigenstates of VT. By choosing a proper phase, 
the eigenvalue of VT can be set to unity: 

VT\E±) = V\E±Y = \E±). (3.22) 

One may think that our definition of CPT-inner product is slightly different from the 
one in the literature. In Refs. j2|, Isj the CPT-inner product was defined as 

(V^, 0)bbj-bmw = {CVTm^\(P) = {CVWfl^) = mV^C^l^), (3.23) 

where T denotes matrix transpose. This definition leads to a weight function VFrbj-emw = 
V^C^. Since both V and C are symmetric in this case, the two definitions of the inner 
product are actually the same. 



D. Special case 4: Mostafazadeh Hamiltonian 
In Ref. jol , Mostafazadeh introduced a five-parameter-class of Hamiltonians with the form 

TT _( r -|- t COS0 — issin0 tsin0 -|- i(scos0 — m) \ ('iOA\ 

^Mostafazadeh ' [^t SIH (j) + i{s COS (j) + u) T-t COS + is siu J' > 

This is a special case of our general Hamiltonian in (13. 6 p by the following parameter mapping: 



e — )■ r, 7 — )• 



a/ sin^ (f) + u"^ a/ sin^ + 

^ t COS0 u ,„ 

cost^ —7- , tan . (3.25) 

Vt2 + t sm 

E. Special case 5: Mostafazadeh-Ozgelik Hamiltonian 



In Ref. [lOj, Mostafazadeh and Ozgelik constructed a six-parameter-class of Hamiltonians 
with a very elegant form 

u -n I z?/^ ^os® e"'*sine\ , . 

HMO = qt,^2 + Ei^^,^^.^^ cose J' (3-26) 

where q and are real and O and $ are complex. In principle, this Hamiltonian is equivalent 
to our Hamiltonian in (13. 6p by the following parameter mapping, 

e q, ±^7^ - /i^ - -> iJ, 

^ ^ 7 sin 6^ + iu cos 6* — z/ r,-„ r,-^ 

7cos^-i/isin0 ^ EcosO, ^ e^^'^ ^ e^'*. 3.27 

7 sm fc* + i/x cos + u 

However, at the degenerate point, E = in (13.260 or 7^ = /i^ + z/^ in (13.60 . these two 
parametrizations are no longer equivalent. In this special point, i^MO is proportional to the 
identity matrix but H in (13. 6p is not. 
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IV. 3x3 CASE 

In this section, we reveal the general form of "PT-symmetric 3x3 matrix Hamiltonians. In 
this case, we use the Gell-Mann matrices, which are the generalization of the Pauli matrices 
in 3 X 3. Any 3x3 matrix can be written as a linear combination of the identity matrix 
and the eight Gell-Mann matrices. If a matrix is Hermitian, the expansion coefficients are 
all real. 

There are two types of solutions for the V operator as well. The first type is the trivial 
solutions Vq = ±1. The second type is the non-trivial solutions. If we expand the non-trivial 
solutions of V as 



P3x3 = ± Pol + J]P.A, , (4.1) 



i=l 

where Aj are Gell-Mann matrices for i = 1, ■ ■ ■ ,8. We find that 

^0 = ^ (4.2) 

and that coefficients Pi depend on four independent parameters. We thus choose the 
parametrization as 

P4 = sin 2x sin 9 cos ip, 
P5 = sin 2x sin 6 sin 
Pq = sin 2x cos 9 cos p, 

P7 = sin 2% cos 6' sin p. (4-3) 

The other four components. Pi, P2, P3, and Pg depend on the sign of cos2x. For the case 
of cos2x > 0, we have 

Pi = — sin^ X sin 26* cos(p — (p), 
P2 = sin^ X sin 26' sin(p — (f), 
P3 = sin^ X cos 29, 

P8 = ^(l + 3cos2x). (4.4) 

Plugging Pq and Pi into (14. ip . we get the non-trivial four-parameter solutions of the V 
operator with the form 

cos 2x sin^ 9 + cos^ 9 - sin^ x sin 29 e'^^"'^^ sin 2% sin 9 e"'*^ 
P3X3 = ± I - sin^ X sin 2^ e-^(^-^) cos 2x cos^ + sin^ ^ sin 2x cos ^ e"''' |. (4.5) 
sin 2x sin 9 e"^ sin 2% cos 9 e''' — cos 2% 

For the case of cos 2% < 0, the correct results are obtained by replacing cos 2% by — cos 2% 
in (113!) and (S3]). 

Just like in the 2x2 case, PT-symmetric Hamiltonians can be found by solving an 
eigenvalue problem. If we use the ansatz of the Hamiltonian as 

8 
i=l 
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the VT symmetry leads to the conditions of e = e* and that ctj satisfy the eigenvalue 
equation 



^Mfc,ai = a*, (4.7) 



1=1 



where 

Mfc, = P^5ki + 2Po Pj^''' + + H ^^-^"^ (^'''^'"' + r7""') (4.8) 

j jml 

with c?*-''^ and /'■'^ being symmetric and antisymmetric structure constants of the SU(3) 
group. 

Once again, the real part of ctj forms eigenvectors of M with eigenvalue +1, and the imag- 
inary part of ai forms eigenvectors with eigenvalue —1. There are always four eigenvectors 
with eigenvalue +1 and four eigenvectors with eigenvalue —1. 

It is straightforward to show that the vector Pi with components defined in fl4.3l) and 
fl4.4p is an eigenvector of M with eigenvalue +1. All other eigenvectors can be constructed 
by the derivatives of Pi. The set of all four first order derivatives, {d^Pi,d0Pi,dpPi,d^Pi}, 
forms a subspace with eigenvalue —1. The remaining three eigenvectors with eigenvalue +1 
can be constructed from the second derivatives. Below is a choice of orthonormal set of 
eigenvectors with eigenvalue +1, 



4^ = —P^ 

2 



73 



— — sin^ X sin 26^ cos(p — </?), sin^ x sin 26* sin(p — </)), sin^ x cos 26', 



2 



1 ~|~ 3 cos 2^ 

sin 2x sin 9 cos y?, sin 2% sin 9 sin sin 2% sin 9 cos p, sin 2% sin 9 sin p. 



2^3 



^ i^2p 3 

« 2-^ 2 



1 3 + cos2x . ( ^ 3 + COS 2^ . w ^ 3 + cos 2x 

sm 2t^ cos(p — Lp) , sm 29 sm(p — j , cos 29, 



2 \ 2 " 2 " 2 

— sin 2x sin 9 cos — sin 2% sin 6' sin — sin 2% cos 9 cos p, — sin 2% cos 6' sin p 

V3 . 2 

^sm X 

= ^— cos X cos 26'cos(p — v^), cos X cos 26'sin(p — (/?), — cos X sin 26*, 

— sin X cos 6' cos ip, — sin x cos 9 sin sin x sin 9 cos p, sin x sin ^ sin p, 
A^'^^ = ^cos X sin(p — ip) , cos x cos(p — (/?), 0, — sin x cos 6* sin yj, sin x cos cos ip, 

sin X sin 6^ sin p, — sin x sin 9 cos p, j . (4.9) 
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Likewise, a set of eigenvectors with eigenvalue —1 can be chosen as 

= - ^— sin 2% sin 26 cos(p — ^p) , sin 2% sin 29 sin(p — Lp) , sin 2% cos 29, 2 cos 2% sin 9 cos 99, 
2 cos 2% sin 6* sin y9, 2 cos 2% cos 6* cos p, 2 cos 2% cos 6* sin p, — a/S sin 2% 

B?" = ^deP. 
2 sm X 

- sin X cos 26^ cos(p — ip) , sin x cos 29 sin(p — </)), — sin x sin 29, 
cos X cos 6^ cos ip, cos x cos 9 sin — cos x sin 6^ cos p, — cos x sin 6' sin p, 

^(3) _ dcf,Pi + 9<^Pi 



sin 2x 

0, 0, 0, — sin 9 sin (p>, sin 6* cos (p>, — cos 6* sin p, cos 6^ cos p, 
i?!^'' = ^sin X sin(p — sin X cos(p — 0, cos X cos 6* sin 

— cos X cos 6^ cos p>, — cos x sin 9 sin p, cos x sin 6^ cos p, ) . (4-10) 



With the solution of ai = Ai + iBi, we can construct the general 3x3 PT-symmetric 
Hamiltonian by plugging the above results into the ansatz in (14. 6p : 



i=l 



(4.11) 

This construction has thirteen real parameters, four in V, four 7's, four p's, and one e. 
Any 3x3 Hermitian Hamiltonian can be considered as a special case with all p's vanishing, 
which also has the correct number of parameters: nine. Unlike the 2x2 case, the general 
PT-symmetric matrix Hamiltonian in (14. lip does not present all 3 x 3 matrices with all 
real eigenvalues. This can be seen by a simple parameter counting. A 3 x 3 matrix with 
all real eigenvalues should have fifteen real parameters but -ffsxs in (14. lip only has thirteen 
parameters. 



V. CONCLUSION 

In this paper, we find the general V operator and construct the general PT-symmetric 
matrix Hamiltonians in 2 x 2 and 3x3. In both cases, VT symmetry can be considered as 
a generalization of Hermiticity. We conjecture that this statement is true in general. 

We convert the searching for a general PT-symmetric Hamiltonian problem to an eigen- 
value problem. This method also applies to higher dimensions. For example, the definition 
for the matrix M^j in (14.81) can be easily generalized to N dimensions by replacing 3 by N. 
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